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ABSTRACT 

The cluster correlation function and its richness dependence are determined 
from 1108 clusters of galaxies - the largest sample of clusters studied so far 

- found in 379 deg^ of Sloan Digital Sky Survey early data. The results are 
compared with previous samples of optically and X-ray selected clusters. The 
richness- dependent correlation function increases monotonically from an average 
correlation scale of ~ 12 Mpc for poor clusters to ~25 Mpc for the richer, 
more massive clusters with a mean separation of ~90 Mpc. X-ray selected 
clusters suggest slightly stronger correlations than optically selected clusters (~ 
2-0"). The results are compared with large-scale cosmological simulations. The 
observed richness-dependent cluster correlation function is well represented by 
the standard flat LCDM model (fim— 0.3, h —0.7), and is inconsistent with the 
considerably weaker correlations predicted by fim= 1 niodels. An analytic rela- 
tion for the correlation scale versus cluster mean separation, tq — d, that best 
describes the observations and the LCDM prediction is tq ~ 2.Q^/d (for d ~ 20 

- 90 Mpc). Data from the complete Sloan Digital Sky Survey, when avail- 
able, will greatly enhance the accuracy of the results and allow a more precise 
determination of cosmological parameters. 

Subject headings: cosmology:observations-cosmology:theory-cosmological parameters- 
dark matter-galaxiesxlusters: general-large-scale structure of universe 
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1. Introduction 

The spatial correlation function of clusters of galaxies and its richness dependence pro- 
vide powerful tests of cosmological models: both the amplitude of the correlation function 
and its dependence on cluster mass/richness are determined by the underlying cosmology. 
It has long been shown that clusters are more strongly correlated in space than galaxies, by 
an order-of-magnitude: the typical galaxy correlation scale, ~ 5 h"^ Mpc, increases to ~ 
20 - 25 h'^ Mpc for the richest clusters (Bahcall & Soneira 1983; Klypin & Kopylov 1983; 
see also Bahcall 1988; Huchra et al. 1990; Postman, Huchra, & Geller 1992; Bahcall & West 
1992; Peacock & West 1992; Dalton et al. 1994; Croft et al. 1997; Abadi, Lambas, & Muriel 
1998; Lee & Park 1999; Borgani, Plionis, & Kolokotronis 1999; Collins et al. 2000; Gonza- 
lez, Zaritsky, & Wechsler 2002; and references therein). Bahcall & Soneira (1983) showed 
that the cluster correlation function is richness dependent: the correlation strength increases 
with cluster richness, or mass. Many observations have since confirmed these results (refer- 
ences above), and theory has nicely explained them (Kaiser 1984; Bahcall & Cen 1992; Mo 
& White 1996; Governato et al. 1999; Colberg et al. 2000; Moscardini et al. 2000; Sheth, 
Mo, & Tormen 2001). But the uncertainties in the observed cluster correlation function as 
manifested by the scatter among different measurements remained large. 

In this paper we use the largest sample of clusters yet investigated, 1108 clusters selected 
from 379 deg^ of early Sloan Digital Sky Survey data (see the SDSS cluster catalog: Bahcall 
et al. 2003b), to determine the cluster correlation function. This large, complete sample of 
objectively selected clusters, ranging from poor to moderately rich systems in the redshift 
range z = 0.1 - 0.3, allows a new determination of the cluster correlation function and its 
richness dependence. We compare the SDSS cluster correlation function with results of pre- 
vious optically and X-ray selected clusters (§3). We use large-scale cosmological simulations 
to compare the observational results with cosmological models (§4). The data are consistent 
with predictions from the standard fiat LCDM model (Q^~0.3, h ~0.7), which best fits nu- 
merous other observations (e.g., Bahcall, Ostriker, Perlmutter, & Steinhardt 1999; Bennett 
et al. 2003; Spergel et al. 2003). 

2. SDSS Cluster Selection 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) will provide a comprehensive 
digital imaging survey of 10^ deg^ of the North Galactic Gap (and a smaller, deeper area 
in the South) in five bands {u, g, r, i, z) to a limiting magnitude of r <23, followed by 
a spectroscopic multi-fiber survey of the brightest one million galaxies, to r <17.7, with a 
median redshift of z ~0.1 (Pukugita et al. 1996; Gunn et al. 1998; Lupton et al. 2001; Hogg 
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et al. 2001; Strauss, et al. 2002). For more details of the SDSS see Smith et al. (2002); 
Stoughton et al. (2002); and Pier et al. (2003). 

Cluster selection was performed on 379 deg^ of SDSS commissioning data, covering the 
area q;(2000) = 355° to 56° and 145.3° to 236.0° at 5(2000)= -1.25° to 1.25° (runs 94/125 
and 752/756). The clusters studied here were selected from these imaging data using a color- 
magnitude maximum-likelihood Brightest Cluster Galaxy method (maxBCG; Annis et al. 
2003). The clusters are described in the SDSS cluster catalog of Bahcall et al. (2003b). The 
maxBCG method selects clusters based on the well-known color-luminosity relation of the 
brightest cluster galaxy (BCG) and the E/SO red ridgeline. The method provides a cluster 
richness estimate, Ng^i (the number of E/SO galaxies within 1 h^^ Mpc of the BCG that are 
fainter than the BCG and brighter than Mj(hm) — -20.25), and a cluster redshift estimate 
that maximizes the cluster likehhood (with l-cr uncertainty of Uz — 0.014 for Ngai>lQ and 
Gz — 0.01 for Ngai>20 clusters). We use all maxBCG clusters in the estimated redshift 
range Zest = 0.1 - 0.3 that are above a richness threshold of Ngai> 10 (corresponding to 
velocity dispersion >350 km s^^); the sample contains 1108 clusters. The selection function 
and false-positive detection rate for these clusters have been estimated from simulations and 
from visual inspection to be < 10% (Bahcall et al. 2003b). 



3. The Cluster Correlation Function 

The two-point spatial correlation function is determined by comparing the observed 
distribution of cluster pairs as a function of pair separation with the distribution of random 
catalogs in the same volume. The correlation function is estimated from the relation ^cc{f) = 
Fdd{''^) / FRR{r) — 1, where F£,£){r) and Fji^i^r) are the frequencies of cluster-cluster pairs as 
a function of pair separation r in the data and in random catalogs, respectively. The random 
catalogs contain ~ 10^ times the number of clusters in each data sample; the clusters are 
randomly positioned on the sky within the surveyed area. The redshifts of the random 
clusters follow the redshifts of the observed clusters in order to minimize possible selection 
effects with redshift. Comoving coordinates in a flat LCDM cosmology with ^Im— 0.3 and a 
Hubble constant of Hq — 100 h km s~^ Mpc~^ are used throughtout. 

The uncertainty in the estimated cluster redshifts {az — 0.01 for Ngai>20 clusters and 
az — 0.014 for Ngai>10 to >15 clusters; §2) causes a small smearing effect in the cluster 
correlations. We use Monte Carlo simulations to correct for this effect. We use simulations 
with a realistic cluster distribution with redshift and richness, convolve the clusters with the 
observed Gaussian scatter in redshift as given above, and determine the new convolved cluster 
correlation function. As expected, the redshift uncertainty causes a slight weakening of the 
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true correlation function, especially at small separations, due to the smearing effect of the 
redshift uncertainty. We determine the correction factor for this effect as a function of scale r 
from 10^ Monte Carlo simulations for each sample. The correction factor (typically < 20%) 
is then applied to the correlation function. An additional small correction factor due to 
false-positive detections is also determined from Monte Carlo simulations using the estimated 
false-positive detection rate of 10%±5% for Ngai>lO clusters, 5%±5% for Ngai>l3, and < 5% 
for the richest clusters with Ngai>15. The correlation function uncertainties are determined 
from the Monte Carlo simulations. Each simulation contains the same number of clusters as 
the relevant data sample. The final uncertainties include the statistical uncertainties and the 
uncertainties due to the small correction factors in the redshift and false-positive corrections. 

The correlation function is determined for clusters with richness thresholds of Ngai> 10, 
>13, >15, and >20. The space densities of these clusters, corrected for selection function 
and redshift uncertainty (Bahcall et al. 2003a), are 5.3 x 10-^ 2.2 x 10"^ 1.4 x 10-^ and 
0.5 X 10~^ Mpc~^ {zest = 0.1 - 0.3). The correlation function of the four samples are 
presented in Figure 1 and Table 1. The best-fit power-law relation, ^(r) = i^J^^- derived 
for r < 50 Mpc, is shown for each sample. The power-law slope 7 has been treated both 
as a free parameter and as a fixed value (7 = 2). The difference in the correlation scale Tq 
for these different slopes is small (< 2%), well within the measured uncertainty. 

The richness dependence of the cluster correlation function is shown in Figure 2; it 
is presented as the dependence of the correlation scale Tq on the cluster mean separation d 
(Bahcall & Soneira 1983; Szalay & Schramm 1985; Bahcall 1988; Croft et al. 1997; Govcrnato 
et al. 1999; Colberg et al. 2000). Samples with intrinsically larger mean separations corre- 
spond to lower intrinsic cluster abundances {rid = d^^) and thus to higher cluster richness 
and mass (for complete samples). We compare our results with those of previous optically 
and X-ray selected cluster samples (Figure 2). These include the correlation function of 
Abell clusters (Bahcall & Soneira 1983; Peacock & West 1992; Richness class >1; Richness 
= clusters are an incomplete sample and should not be included); APM clusters (Croft et 
al. 1997); Edinburgh-Durham clusters (EDCC; Nichol et al. 1992); Las-Campanas Distant 
Cluster Survey (LCDCS; Gonzalez, Zaritsky, & Wechsler 2002); galaxy groups (2dF; Zandi- 
varez, Merchan, & Padilla 2003); and X-ray selected clusters (REFLEX: Colhns et al. 2000; 
and XBACS: Abadi, Lambas, & Muriel 1998; Lee & Park 1999). A summary of the results 
is presented in Table 1. For proper comparison of different samples, we will use the same 
set of standard parameters in the relative tq - d plot: redshift z ~0, correlation power-law 
slope 7 = 2, and all scales are in comoving units in the LCDM cosmology. We discuss these 
below. 

Most of the cluster samples are at small redshifts, 2; < 0.1 (Table 1). The only exceptions 
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are the SDSS clusters ~ 0.1 - 0.3), and the LCDCS {z ~ 0.35 - 0.575). To convert the 
results to 2; ~ we use large-scale cosmological simulations of an LCDM model and determine 
the cluster correlation function and the Tq - d relation at different redshifts. Details of the 
simulations and cluster selection are given in Bode ct al. (2001) (sec also §4). The correlation 
function is determined following the same method used for the data. We find that while 
both To and d increase with redshift for the same mass clusters, as expected, there is no 
significant change (< 3%) in the - d relation as the redshift changes from 2; = to ~ 0.5 
(for d ~ 20 — 90 h"^ Mpc). In Figure 2 we plot the individual parameters Tq and d at the 
sample's measured redshift as listed in Table 1; the relative ro- d relation remains essentially 
unchanged to 2; ~ 0. 

Most of the cluster correlation functions (Table 1) have a power- law slope in the range 
of ~ 2 ±0.2. The APM highest richness subsamples report steeper slopes (3.2, 2.8, 2.3); they 
also have the smallest number of clusters (17, 29, 58). The correlation scale ro is inversely 
correlated with the power-law slope; a steeper slope typically yields a smaller correlation 
scale. We use the APM best fit foi" at 7 = 2 (Croft ct al. 1997) for these richest 
subsamples. Using cosmological simulations, we investigate the dependence of tq on the 
slope within the more typical observed range of 2 ± 0.2. For the current range of mean 
separations d we find only a small change in tq (^5%) when the slope changes within this 
observed range. The X-ray cluster sample XBACS yields similar correlation scales for slopes 
ranging from ~ 1.8 to 2.5 (Abadi, Lambas, & Muriel 1998 and Lee & Park 1999). Similarly, 
the SDSS correlation scales are essentially the same when using a free slope fit (typically 
1.7 to 2.1) or a fixed slope of 2. Since most of the observations are reported with a slope 
of 2, we adopt the latter as the standard slope for the results presented in Figure 2. The 
only correction applied is to the three highest richness APM subsamples; these are shown 
both with and without the correction. We also varify using cosmological simulations that 
the LCDCS sample at 2; ~ 0.35 - 0.575, with a slope of 2.15, has anvQ - d relation consistent 
with the standard set of parameters used in Figure 2 (2; ~ 0, 7 ~ 2). 

Finally, we convert all scales (rg and d from Table 1) to the same r2m,=0.3 cosmology 
(LCDM). The effect of the cosmology on the observed Tq - d relation is small (< 3%), partly 
due to the small redshifts, where the effect is small, but also because the cosmology affects 
both To and d in the same way, thus minimizing the relative change in the Tq - d relation. 

A comparison of all the results, including the minor corrections discussed above, is 
shown in Figure 2. Figure 2a presents both optically and X-ray selected clusters; Figure 
2b includes only the optical samples. The richness-dependence of the cluster correlation 
function is apparent in Figure 2. The X-ray clusters suggest somewhat stronger correlations 
than the optical clusters, at a ~2-(j level. Improved optical and X-ray samples should reduce 
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the scatter and help address this important comparison. 

4. Comparison with Cosmological Simulations 

We compare the results with large-scale cosmological simulations of a standard LCDM 
model 0.3, h = 0.67, (jg = 0.9), and a tilted SCDM model, TSCDM 1, h ^ 0.5, 

n — 0.625, (78 = 0.5). The TPM high-resolution large- volume simulations (Bode et al. 2001) 
used 1.34 X 10^ particles with an individual particle mass of 6.2 x 10^^h~^ M©; the periodic 
box size is 1000 h'^ Mpc for LCDM and 669 h'^ Mpc for TSCDM. The simulated clusters 
are ordered by their abundance based on cluster mass within 1.5 Mpc. The results of 
the cosmological simulations for the Tq - d relation of z = clusters are presented by the two 
bands in Figure 2 {1-a range). A correlation function slope of 2 was used in the analysis. 
We also show the simulations results of Colberg et al. (2000) for LCDM, and Governato et 
al. (1999) for a standard untilted SCDM (fi^= 1, /i = 0.5, n = 1, (Xg = 0.7). The agreement 
among the simulations is excellent. As expected, the untilted SCDM model yields smaller 
ro's than the strongly tilted model; LCDM yields the strongest correlations. 

We determine an analytic relation that approximates the observed and the LCDM ro - 
d relation: ro ~ 2.6Vd (for 20 < d < 90; all scales are in Mpc). The observed richness- 
dependent cluster correlation function agrees well with the standard LCDM model. The 
correlation scales, and the - d relation, increase as Jl^/i decreases and the spectrum shifts 
to larger scales. The flm= 1 models yield considerably weaker correlations than observed. 
This fact has of course been demonstrated earlier; in fact, the strength of the cluster cor- 
relation function and its richness dependence were among the first indications against the 
standard n„,= 1 SCDM model (Bahcall & West 1992; Croft et al. 1997; Borgani, Phonis, & 
Kolokotronis 1999; Governato et al. 1999; Colberg et al. 2000; and references therein). 

The scatter in the observed Tq - d relation among different samples is still large, es- 
pecially when both the optical and X-ray selected samples are included. A high-precision 
determination of the cosmological parameters cannot therefore be achieved at this point. 

5. Conclusions 

We determine the cluster correlation function and its richness dependence using 1108 
clusters of galaxies found in 379 deg^ of early SDSS data. The cluster correlation function 
shows a clear richness dependence, with increasing correlation strength with cluster rich- 
ness/mass. The results are combined with previous samples of optical and X-ray clusters. 
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and compared with cosmological simulations. We find that the richness-dependent cluster 
correlation function is consistent with predictions from the standard flat LCDM model (fim= 
0.3, h = 0.7), and, as expected, inconsistent with the weaker correlations predicted by flm= 
1 models. We derive an analytic relation for the correlation scale versus cluster mean separa- 
tion relation that best describes the observations and the LCDM expectations: tq ~ 2.6Vd. 
X-ray selected clusters suggest somewhat stronger correlations than the optically selected 
clusters, at a ~2-(T level. 
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Fig. 1. — The SDSS cluster correlation function for four richness thresholds {Nyai>10, >13, 
>15, >20). Best-fit functions with slope 2 and correlation-scale Tq are shown by the dashed 
lines (1-cr uncertainties). 
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Fig. 2. — Correlation-scale Tq versus mean cluster separation d for all samples (Fig. 2a) and 
for optical samples (Fig. 2b) {1-a uncertainties). A slope 7= 2 and LCDM comoving scales 
are used (due to conversion to the standard LCDM cosmology, some values differ slightly 
from Table 1; see §3). Cosmological simulations are presented by the two bands (LCDM 
and Tilted-SCDM). Previous simulations of LCDM (dot-dash) and untilted SCDM (dotted 
curve) are also shown (§4). The solid curve represents Tq = 2.6y/d (§4). 
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Table 1. The Cluster Correlation Function 



Sample* 




z 


7 


ro 


d 










(Mpc) 


(Mpc) 


SDSS 














1108 


0.1-0.3 


2 


12.7 ± O.O 


26.6 


Ngai >13 


472 


0.1-0.3 


2 


15.1 ± 0.9 


35.6 


Ngai >15 


300 


0.1-0.3 


2 


17.3 ± 1.3 


41.5 


Ngai >20 


110 


0.1-0.3 


2 


21.2 ± 2.8 


58.1 


Abell (1,2) 












Rich >1 


195 


<0.08 


2 


21.1 it 1.3 


52 


APM (3) 














OD4 


n 1 \ 

^<U.i> 


O 1 


1 /I o -1-0.4 
14.2 ±0 g 


oU 


R>70 


114 


<0.1> 


2.1 


16.6 ± 1.3 


48 


R>80 


110 


<0.1> 


1.7 


18.4 ±2;2 


57 


R>90 


58 


<0.1> 


2.3 


22.2 ± 2.8 


69 








[ 2 ] 


[23.0 di 2.9]^ 




rVj^ ±uu 








±(3.4: ^ 4.(3 


7Q 








r o 1 
[ 2 ] 


[22.0 ± D.OJ 






1 V 


<U.1> 


O.Z 


1 Q _1_ IT Q 


OO 








a\ 


i27.0 ± ().7]'' 




EDCC (4) 


79 


^0.13 


2 


16.2 ± 2.3 


41 




178 


0.35-0.475 


2.15 


14.7 ±2 2 


38.4 




158 


0.35-0.525 


2.15 


17.2 ±\% 


46.3 




115 


0.35-0.575 


2.15 


on fi 1 4.6 

20-9 ±6.6 


58.1 


REFLEX (6) 












Lx >0.08 


39 


£0.05 


2 


24.8 ±2.5 


48 


Lx >0.18 


84 


<0.075 


2 


25.8 ±1:9 


61 


L,Y >0.3 


108 


<0.10 


2 


31.3 ±2;? 


72 


Lx >()••') 


101 


<0.125 


2 


25.8 =:|-:| 


88 


XBACS (7,8) 












Ljf >0.24 


49 


£0.07 


1.8-2.5 


25.7 ±3.7 


66 


Lx >0.48 


67 


£0.09 


1.8-2.5 


25.2 ±4.1 


82 


Lx >0.65 


59 


<0.11 


1.6-2.2 


30.3 ±i:i 


98 


Groups (9) 












>5cl3 


920 


<0.12> 


2 


11.7±0.6 


20.9 


M„ >lel4 


540 


<0.13> 


2 


13.4 ±0.9 


28.9 



'^Sample (with reference), and subsample threshold in richness, X- 
ray luminosity (10'*'* erg s^-*^), or M^i^ (Mq). References: l.Bahcall 
k. Soneira 1983; 2.Peacock & West 1992; 3.Croft ct al. 1997; [larger ro's 
are obtained for APM by Lee & Park 1999]; 4.Nichol et al. 1992; S.Gon- 
zalez, Zaritsky, & Wechsler 2002; 6.Collins et al. 2000; 7.Lee & Park 
1999; see also 8.Abadi, Lambas, & Muriel 1998; 9.Zandivarez, Merchan, 
& Padilla 2003. The SDSS, LCDCS, and Groups use LCDM cosmology 
for their ro and d; all others use f2rn=l. AH scales are for h = 1. 

''Correlation-scale ro using a slope of 2 (see §3) 



